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ifflSTRACT 

{Pheoretical perfariaance of mltistage axial-flow turbomachinos 
is determined using the streamline curvature method. She mathematical 

7 

model is based on that developed by Moorthy , !Ehe flow field equations 
are solved numerically on a computer. Input to the program consists 
of annulus geometry, design fluid angles, a poly tropic coefficient 
to account for losses at each operating point and inlet conditions 
to the turbomachine. 

Ehe radial equilibrium equation (r. e. e. ) along with continiiity 
equation and auxiliary equations is solved in the inter blade row 
gaps using finite difference technique% The output consists of 
radial distribution of flow parameters, overall pressure ratio and 
adiabatic efficiency as a function of mass flow and rotational speed. 
These results are compared with experimental resxrlts, whenever 


available^ 
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NOMEUCLiTUEE 


A 

■^1 


B 

G 

C 

C 

■P 

D 

e 

f 

H 

h 


Defined in Eq;,- H-'? 

Constants used in estination of deviation angles at off' 
design conditions 

Defined in Eq. II -9 

Defined in Eq, II -.10 

Absolute Velocity 

Specific heat at constant pressure 

Equivalent diffusion factor 

Fraction of mss flow through each streantube 
Stagnation, enthalpy 
Static enthalpy 


M 


M 


abs 


M 


rel 




n 

V 

<1 

r 

s 

T 

t 

U 


Mober of calculating stations 
Meridional Mach number 
Absolute Mach Mimber 
Relative Mach number 
RUmber of streamlines 
Polytrop^c ooeffioient 
static pressure 

Any flow quantity defined in Eq. II-7 

Radius 

Entropy 

Stagnation tempemture 
Static temperature 
Rotor velocity 




Mass flow through the turhomachine 


w 

Relative velocity 


X 

Axial distance 


a 

Ahsolute fluid angle measured with respect to 
direction 

axial 

e 

Relative fluid angle measured with respect to 
direction 

axial 

A 

Defined in Eq, 11-17 


fi 

Deviation angle 


n 

Stage efficiency 


X 

Meridional streamline slope 


Xt 

Meridional streamline curvature 

- 


Relaxation factory on axial velocity 


^2 

Relaxation factor on streamline shift 


(0 

Angular velocity of rotor 


p 

Density 

' ' ' ' . ' ' ' ■ i 

Y 

Ratio of specific heats 

■ ' 't' 

Z 

Defined in Eq, 11-18 

■ " 'I 

3x113303: ip ts 

■ ■■ 

i 

■ ■ ■ ■ ■ ' I 

o 

Inlet or Stagnation conditions 

' i 

I 

o 

Compressor 

■ - . ■■ ■ '.i 


Denotes calculating station 

I 

in 

Inlet to a hlade row 


d 

Denotes streamline numhec 


r 

Component along radial direction 

!, 

t 

turhine 




X Component along axial direction 

e ■ Component along tangential direction 

Superscripts : 

• Relative to rotor 

* Design or minimum loss values. 



CHAP‘IER I 


IlgROIiIJCTIQI? 

I.A. GEEBML BACKGROUND 

Before finalising the design of a tarhomachine, a designer has 
to consider various annulus geometries and hlade geometries. He has to 
consider design performance as well as performance at off-design 
conditions. Often a machine performing well at design condition say 
give quite unsatisfactory performance at off-design conditions. Hence 
thfonetical methods of performance prediction of turboiiia.chine heoome 
neoeesary and by employing such methods various designs can be studied 
for their performance at all operating conditions. A design can thus 
be finalised with some certainty before actual manufacture and testing 
with their problems of high cost and time delays. 

Theoretical approach to the solution of aerodynamio problems 

in the design or perfomance prediction of fluid flow in axial turbo- 

machines presents formidable problems. It involves the solution of the 

usual non-»linear differential equations of fluid flow' in most complicated 

geometries. An exact solution of the flow field, governed hj^- the 

equations of three-dimensional unsteady, nordsentrcpio viscous 

flow through the complex of rotating and stationary hladcs and 

associated ducts is virtually impossihle and hence calculation 

procedures must he approximate. Most effort has been concentrated 

on the formulation of meaningful two-dimonsxc ml models oi Xiie 
: / V I 

flow field viz , 1 ) the hub-to-casing or meridional flow- 

model .• e.g. the S2 s-urface of Wu"* and 2) the blade-to-blade or cascade 



3 


1 

mod^ e*g. the Si surfece of Wu • bsisis of each laodel is a Judicious 
mixture of theoretical technique experimeiibal correlations end experience* 
In the present insrestigation only sttb3y of hub-to-tip analysis or meridional 
flow model is undertaken* The fluid angles and fluid turning losses, which 
ere necessary for the performance prediction of turbdmechines ere usually 
obtained from the cascade flow model. In the present case design fluid, 
angles are fed in as input and the losses are accounted for by introduction 
of 8 polybropic eKponart of expansion or ccnpression d^ ending on whether 
the machine is a turbine or a conpressor* 

In the earlier studies on turbomachines, th'^' flow was assumed 
to be cylindrical and onHy velocity conponents in axial and tangential 
directions were considered* The radial conponent of velocity was neglected* 
This was quite reasonable with the earlier type of turbomachines where the 
blnde height wes usually ssnaH, the hubMnp ratio was high and the stream- 
line were usually parallel to the axis* Bub xdth the advent of modern Jet 
engine compressors add large steam and gas turbines, the hub-tip ratio has 
become low and the blade height is significant* Hence the flow tends 
towards conical rather than cylindrical and iiie radial ccraponenb of velocit: 
becomes significant* SfareaioLine d-ope and curvature, therefore, play a 
dominanb role in the flow field analysis* The ^ape of a streamline has a 
great effect on the rediel pressure gradient and the radial distribution 
of other flow properties* The radi^ dii^xLacement of the fluid introduces 
curvjjibure into the streamlines and this introduces additional radial 
pressure gradient, which is stpeiinposed on the pressure gradient due to 
tangential motion. The changes in the radial pressure component can be 
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qm-te large in the case of long blades, i»e* the low pressure end of a 
turbomachine, even though the actual radial displacement of the stream- 
line may be small. The magnitude of radial di^la cement of streamline 

increases considerably with the Ifeeh nunber. 

^ctca-nt 

In order to take into^the effects of i^dial di^lacemenb of stream-- 
lines, the streamlines through the stage must be determined as accurately 
as possible for perfoimanee calciiLations. The slope of the streamlines 
allows one to determine the radial velocity conponent and the cuivatura 
enables the radial pressure gradient to be calcuLaited. Thus, the slope 
and curvature of sbreanlinB becomes -•important in the performance calcu- 
lation of low hub-tip ratio turboma chine and high Ifech nunber flows* 

In the past the methods that have been employed for calculating 
the flow field in axial turboma chines were usually based on relaxation 
techniques, strsamLime—curvature (s.C*) methods, actuator disk theoirf-ss, 
siEple radial equilibrium theories or matrix methods, Becenfcly with the 
advent of fast digital computers and improved conputational techniques, 
streamline curvature and matrix through flow methods have been favoured* 

Tn the present investigation the streamline curvature db thod is 
used to predict the performance of an axial turbomachine when amulus 
geomietry and blade geometry ard a polybropic exponent to account for 
lossesi turbomachine inlet conditions such as stagnation 

pressure and t®iiperatU3^, which are assumed t(^ constant raaiaJij^ are 
given. The radially distributed temper eture rise, pressure ratio and 
efficiency are subsequently determined as a function of m^ss flow arid 
rotor ^eed* 
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In the S»C» method 0 curve is fitted through the points of 
intersection of slaream surfsces wii^ the radial lines and the properties 
of the curve, i«e» slope and curvature are taken to represent the sfcraanf 
line properties* These properties ore used to calculate radisl component 
of velocity end radial pressure gradient# 

Before going into the details of the present study, a reriaw of 
iuportant work dons previously is presented in the neoct section# 


I.B. RETIEtJ OF EREg-IODS WCRK; 

Some of the earliest work dene on the performance prediction of 

2 3 

axial turbomachines was reported by Seoxvy and inderson and by S-^isnr 

whose work we.s based on the stream filament theory of Ho!U)toquist and 

Hannie'^* Both investigations were based on a finite-differer»e solution 

of the nonisentropic radial equilibrivm and continui-ty equations at 

stations between blade rows. The steady, axii^nmnetric model of the flow 

4,5,6,7,S,9,10 

Used has been described thoroughly by later workers * 111 

each method, correlations of experimaibal data were developed and used 
for predicting the radial distribution of the fluid turning angle 
and of total-pressure loss for each blade. row# Htamarical • solutions for 
single-stage geometa'ies were presented, however, conparison with 
sxpurimental resc3.ts was not good, probably due to inadequate data 
correlations for flow angle and losses# Nevertheless, thejaSSF studios 
denonstarated the feasibility of mmerical solution ^sterns in generating 
both detailed flow passage distributions and overall conpressor 


performaaaoe# 
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A similar approach in dev-eloping a program for conpubstion 

11 

of muLtisfcage sctial-flow oonpressors was described by Jansen and Mjff att • 

12 

Iho steady, asd. symmetric model i^as^gain used as formulated bs’- Ifovak 
and this formulation included iuproved schemes for iterative location 
of the jssi^nmaetric stream surfaces and for computation of local values 
of stream surface slope and curvature. 

Recently 0ravelli3ig and Carmody also have developed an 
analysis program for muLtisbage axisl—flow compressors* The documenta— 
tion of the computer program is incoiplete* In reference 14 Ifeneshlg**^ 
and Shalaan used a flow model similar to Novak and their own data cor- 
relations to predict ccmpressor performance. They discuss in someiirfiat 
more detail than in aiy earlier work, the numerical problems which are 
encountered in flow field solutions* They also examine the ecd-stence and 
uniqueness of the solutions obtained. This useful discussion is st^ple- 
menbed ty papers of snd mUdLnson’^^ who give a good deal of ! 

insight into some of these numerical difficulties* More recently, 

Magumdar ” describes a program developed by Davis and «oiippares it with 

T 12- 

others work* This program is primarily based on the flow model of Novak 
combined with correlations available from earlier studies* 

/m of the preceding work is concsmsd with analysis for ax-ial“ 
flo'vrconpressor geometries* These methods all assume a steady, axi^jr* 
mmetric flow. 4L1 osculations are at stations located in the axiS 
spaces between blade vows* The basic flow modS is, therefore, that 
identified as the blade-element modS and described in reference 10* 
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Tiiisted "burboine chine blades ere asually gyclved by a radial stacking 
of a ssnes of blade sections or blade elements ■which are "breabed 
locally as 2*”D sections to fonn. complete blade* The blade“'eleiQ 3 nb 
model^ hence^ proposes that the blade*T?ow characteristics can be dater~ 
mined if the performance of each blade element is kno'wn and that Ofvej?— 
all machine performance can then be obtained by an integration of the 
performance of each blade row* Fundamental differences in strategy 
exist between the above methods in the numerical solution ■bechniqi^s 
applied, including the means for estimating local vlues of stream 
surface position, slope and curvature* 

A somewhat different means for formulation of the performance 

prediction problem for exial flow turbomachines has been propos ed and 
Used by Mar^ • This method deritas a Poisson type equation for the 
s1a:‘eamf unction* ilso ealeiiLating -planes are taken inside Hade rows, 

20 

by assuming the ^ape of mean stream surface in the blade row* Frost 

comibines this method - the matrix through-flow method - -with the S*C* 

21 

method. Jfcre recently, Ifersh and Bogman have improved the matrix- 
through flow method to precdct performance of various types of turbo- 
machines* iHxis method indudes an improved loss modd based on Horlock*s 
T^k^* The matrix through flow method does not avoid the requirement 

for input of key enpiiical data correlations and in example solutions 
currently available, does not appear to offer a dear inprovemant in 
a^y area of performance calculation except a litile saving in computing 
time* Its disadvantage over the S*C. method is that it requires more 
conputer core storage* 
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Jiddiiiioiisl sx 0 ii|xLas of relslied research on "burbomfichine 
performance prediction j?re conbainsd in references 23 to 25 • i proce~ 
dure »nd eonputer program for axial-flow turtina analysis is described 

o/ 

by Fl'^gg • i33Dther is discussed 1:^ Benaudin and Somm*^. J noreL method 

for avoiding solution oorvergence problems is used by them, which 

25 

should be applicable to other turbomachine cases.Ribaut has outldnead 
a ^stem for a very general analysis of the through-flow field, bub 
unfortunately, the problems in application appear to be substantial* 

In the present insrestigation the flow model is based on the 

7 

formulation developed by M^oithy • This formulation gives a better 
physical picture of the flow and is suitable for rapid digital conpuba- 
-tion. The meridion^ ^ch Number and the streamline slope and curvature 
terms appear explicitly in the radial equilibrium equation, so as to 
provide a better physical appreciation of the effects* of these factors* 
Moreover, when this foimulation is reduced to the two-dimensional case, 
the dynamic equation of fluid flow reduces to the weiLl*dax>wn espression 
for the change in Prandtl-l^er angle acioss a ^fech line. This mbhod 
was employed by Mbortty"^ to study the performance prediction of a 
turbine stage using a high speed digital conputer and the results of 
study wore said to be in very good agreement with experimental values* 

In the next chapters, this method is developed for multistago 
axi?l. turboma chines and developmarib of digital coepuber line up for 
rapid calcuLetion of performance prediction of muLti-sbage turbooBchine 
is preserrbed* 



caapm n 

FCRMOLATICM CP TOE PR.ORTRM 


II.i. PREXJKmiRY CQKSlDERiTIOIB : 

In this chapter the performance cslculation # method for multi-' 

7 

stage turboms chines is developed foUowLng Moorthy • The schemstic dia*" 
grem of meridional flow model is given in Figure 1« The flow is analysed 
at a station i m mediately dovaistream of a blade row by considering the 
known stream conditions, the annulus configurainon and the blade shape* 
This method can be used for both the direct problem, i.e* performance 
prediction of given turbomachine and the inverse problem i.e. determi- 
natLon of cascade characteristics, given other flow properties. The 
present investigation is limited to the direct case. 

Steady, axi^mimatric ccopresable flow is assumed. The celcuLa-r 
ting planes are located in the axial spaces between blade rows. CalcuLa— i 

tions are performed at various radial points at each plane, the radial 

■ . ■ ■ I 

points being selected so that an equal fraction of total mass flow passes 
through each stream tiibe. From the simultaneous solution of -tiie Tadial 
equilibrium and continuity equations, axial velocities and strea mlin e 
radii are determined. Eie effect of streamline ^pe and curvature is 
considered in calculating these quantities, dice axial velocity .and 
Streamline radii are determined, other flow properties are easily deter- 
miaied. The losses due to real fluid effects such as viscosity anal shocks 
are accounted for by the ^ecification of a polytropic oEflgMoaey a cog ffl ^ 
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II. B. iSSUMPIIQNS t 

Since the flow in turboiaechines is quite canplex, some reasonsKLa 
assunptions must be made, ao that s mathematical representation of the 
flow field is tenable. Hence in addition to the previous assumption of 
steady, axissTumetric flow and accounting of losses by introduction of 
polytropic coefficient, the following assimiptions are made. 

1. There is no heat exchange between adjacent stream surfaces. Although 
the tecp era tore difference between root and tip can be quite large in a 
practical machine with long blades, it can be said, that when the fluid 
velocity is large, the adjacent fluid particles do not h^ve sufficient- 
time for any appreciable heat exchange. Further, in the blade-element 
method, the adjacent stream surfaces are assumed to be separated by 
short radial disfc^ces, the actual tenperature difference between them 
is small. Therefore, One may be justified in neglecting heat exchange* 

2» The total or stagnation values of enthalpy, density and angular 
momentum per unit mass of fluid are constant on streamline outside 
any blade row* This does not mean that their values do not change across 
a blade low, i.e. in the passage of the fluid through a row* assumr 

ption means that the gradients of the above properties are subsequently 
zero along the streanUnes* 

3* The boundary lejrer effect is accounted for by the indnsion of 
boundary layer blockage factors, which give the percentage of the anuulus 
area not blocked by boundaiy layers at the hub and the tip. 

4. Interference effects between adjacent blade rows are neglected. The 
flow downstresoa of each blade low is assumed to be mixed completely before 
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entering the next hlade roTr-. 

5 . Tip leakage effects and secondary flo^ effects are neglected 
excent as inclnded in the specification of polytropic coefficient. 
The effects of sweep and dihedral of Hades are neglected. 

6. The specific heat at constant pressure is assumed constant. 


7« Shock losses in transonic compressors are assumed to he 
negligible, 

II. C, DERIVATIOIir 03? HAIIAL EQJJILIBRIUM EQ.IJATIGH s 


In this section, details of the derivation of the radial 

equilihrium eq.uation are given hriefly. This derivation mainly 
7 

follows Moorthy’s formulation* 


The governing equations for steady inviscid axisymmetric 
compressible flow in the absence of body forces and heat transfer 
can be written as follows ; 

Continuity Equation 




(II. 1 ) 


Equation of Motion in Radial direction ; - 

30 3C^ d? . 

'C — ^ + c — — » — « -• — 

X 3x ^ 3r r p 3r 


(II .2) 


Equation of Motion in Axial Pireotign 


3 C 3 C 

„ X ^ X 1 3p 
X 3x r ar P 3x 


(II. 5 ) 



Eaergy Egiiation 



(II.4) 


where C^, C^, are the velocily components, h is the static 
enthalpy and. eq. (II. 4 ) defines H, the stagaation enthalpy. 

CkMabined first and second laws of thermodynamics give 



H 

r 

1 

ar 

ar 

P- ar 


3h 

*• SS5 


3x 

ax 

p ax 


(11.5) 

(11.6) 


where s is the entropy and t, p and p are static values of 
temperature, density and pressure respectively. 

The above equaticaas are now expressed in terms of total 
derivatives along meridional streamlines. Haese streamlines are 
the intersection lines of axisymmetric stream surfaces with tiie 
meridion^ plane containing the r and x directions. Heneefoarth 
streamlines will refer to meridional streamlines. 

Along a streamline, the radius r is an implicit function 
of X. The total derivative of any flow quaitilgr q alaag the 
meridional streamline with reference to x for steady axisymmetric 


flow is given by, 
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or 


dx 


3x 


+ X 


3 

3 r 


(II.9) 


_ dr / dx _ r 
»be re X = — -•?*/« ' — 


dr 

dx 


(II. 10) 


is the streamline slope at the point considered. 

Usisg equation (ll,9j in equations of motion (ll*2 and 
11,3; and combining them after multiplying by 0 and 0 
respectively, we get after some rearrangement, 

r,2 


( 1+ 


dC. 


+ X(C x*<- ) + a 


*x dx 


,2 d / 

dx '• 


log ^ ; = 0 (II. 11, 


d^A> 

where ' a' is the speed of sound and X’ = ■ — ^ is the streamline 

dx*^ 

curvature. The equation of continuity is modified using 
equations (11,9^ and (ll, 10, to 




- + o (-|^ + — ,1 

dx X ' 3 r r 


(II<,12) 


a’liminating — ( log p; from equations (ll,11 and II, -.12/ 
and using the Meridional Mach nuntoe r 


1 , , 2 \ _2 
^ 1+ X ; G 

M , = / — -r — — — 


(II. 13 ) 


equation (ll, 14) is obtained. 


He noe 


dC, 


XM 


dx 


1-M 


1+ X 


rd 


3X X T 

* X* ) + , +' ■ 1: 

3r ^ 


(II. 14' 
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iAs evident from (ll, Inj, for negligible Mach Numbers the 
oircufflfe re ncial velocity and streamline curvature in the 
meridional plane have no effect on axial velocity changes. 

If X and M are both not small, the circumferential ve locity 
and the streamlines curvature becomes important as in the 
case of high subsonic, transonic, supersonic t urb omacbi nes 
with tapered annuli. 

Thus equation gives a batter physical apprscistion 

of effect of Mach number, streamline slope and curvature and 
also indicates that for a singular solution to exist the meri- 
dional Mach number, not the relative Mach number is important. 


Now this equation is used in tbs radial equation of motion 
(ll,2; along with (II.IO; to detain the radial pressure gra- 
dient in a turb omacbi ne , i.e, the radial equilibrium equation. 

Using equation (II, lOj 



dx dx 




X» G 

X 


(II4I5) 


Using this e quation in Squation (ll*2y an.:' substituting for 
dC 

— d? ■■ f roip Squation (H» 1^) obtain 


1 - 

Pp 9^^^^ 


- X C 


1-M 


X M 
1+ X^ 


rC 


2 

X 


- X’J 


3 X 


X 

r 


2 

X’ C - 

X 


2 

Ce 


(II, 16) 
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De f i ni ng 



(II. 17) 


and 


I 



( 1 + 



( 1 + 



(II. 18) 


Equation (II,16> is transf^ngsc; to 


1 3 p 
P 9 r 



1 

2( 1-M^) 



3 r 


+ 


( i+^K 



X' = o (11,19) 


This is the generalised radial equilibrium equation|l^ of 
fcurb omachinis , For turbomachi ne s of cylindrical annuli 
equation# ( 11*19 i becomes simply ” • This equation 

implies that radial velocity is zero at inlet and exit 
from a blade row. For the present investigation neither X nor 
M is negligible, hence A or S are not negligible. 


00 far no statement regarding e ntr opy gs ne ratio n has bee n 
made. In general losses result in a non-uniiorm entropy. How- 
ever it is possible to postulate that the polytropio equation 
P/ ph = constant holds. Using this relation^ Jft equation (11,5/ 
can be writ te n as 


lb = J- 

3r ” Y- 1 n P 9x 


Sow for a compressor 





( 11 , 20 ) 
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and for a turbine 


y n-1 
Y-1 n 


= na 


where n and 'ttie respective stage efficiencies here 

c c 

assumed constant for all stages and for all streamlines. 
Hence 

~ = n — (tuihine) 

3 p t p 3 r 

1^ =! l/x\ — (compressor) 

9r ' 'c p 3r ^ 

Expressing h in terms of H and velocity compaaent the r. e. e. 
for compressors and turbines becomes respectively 


-n 


IH ^ ac f (c") * 1 {„ . -4) oj -5^ 

3r x 2c 9-s 


c 9r 2 


1-M 


2 X 3r 


n sc? o: « 

+ -S — + (l+z) (— - X’ c^) + 
2 3r ^ r 


AO 


(II.21) 


JL IS + 1 itif c2 ^ JL + _JL-} ^ 

ar *2 3r x 2 x 3r 


3C? 
1 6 




« AC 
.2\ X 


( 11 . 22 ) 


II.D. Auxiliary Equations 

In addition to tl^ r.e.e. and ecaatinuity equation a 
number of auxiliary equations are necessaiy to calculate flow 
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in turb omac bines , The continuity equation in the following 
integral form is suitable for the streamline curvature method. 

'li = 2 tt / PC„ r dr (11,23; 

Jim 

^i , 1 ^ nub 


or between two adjacent streamlines the mass flow becomes 




f¥ = 2it / p C r dr 




( 11 , 24 . 


where f is a fraction of mass flow between streamlines (j_i; 
and j, the index i refers to the calculating plane, each of 
which is divided into (iSI-l) stream tubes which are labelled by 
t be i nde x of j . 

The enthalpy rise across a blade row can bs found from 
Ruler's equation 


AH = (lar Qa)i 4 ” ( “r . along streamline j 

W iXp J' J C ^ *** f J 

( 11 . 25 ) 


or 


H. . = H 

1, J 


1-1,0 " (H.a6) 


He ncs te mpe r at ure rise for a perfect gas is 


AT 


( mr C„/. . ■• ( wr Co j . , . 

6 i> J 


IaL 


( 11 . 27 ) 




All the above equatioos are valid for any caloulatlr^ 
plana . Howe ve r the angular ve locity a is non-zero only for a 
rotor, i.e, AH in 11.25 is non-zero only if the plane * i* 
follows a rotor. -Appropriate modifications must also be made 
in the other equations. 

The radial equilibrium equation and continuity equation 
along with the auxiliary equation are solved numerically in 
the inter blade row eap to obtain the axial velocity and 
streamline radius. The me tbod is de scr ±)8 d in detail in the 


next chapter 
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II. E-, gSTII4^TI0.1_ OF ?i.'UI3 AHCxL^o AT 
OFF-paSIGM QOM}ITIO?g } - 

In tbe present investigation, absolute fluid angles 
at the trailing edge station of blade rows at design point 
only are specified. Hence a procedure is necessary to estimate 

3 

fluid angles at other operating conditions. Swan’s method 
of correlation of the deviation angle 6 with respect to 
equivalent diffusion factor is chosen. 

Fig, 3 shows the trend in 6-5* versus Dgq - Dsq* 
which was found from an examination of Double Circular Arc 
(D,C,a,) rotor blade elements. The data were grouped in rela- 
tive i'iach number ranges to show the effect of this parameter. 

No noticeable effect of spanwise location occurred. In general, 
the variation is represented by the following e quation. 

6-6*= ■■ A^ ” ^eq) (ll, 32 j 

where starred quantities represent design values, and A^ = 6.40 
and A^ = 9.45 for P.C.a, rotor blacje elements. 

The constant 0,6 was used because data obtained at relative 
Mach numbers below 0.6 did not indicate any Macb number influ- 
s nee. Similar correlations were cribtained from limited stator 
data. No Mach number effect was applied, because of a lack 
of data on stators in high Mach number range. Thus 

6 - 6 * . A 3 (Etq - 33b*<j) (11.33) 
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where = 2,4 for stator blade slsments. 

Once 6 is de te rtni tjs d , fluid angle at the exit plane 
of a blade row is determined as follows. 


(11.34) 

Due to 1 ask of data, for other blade types, the 


off“design fluid angles are assumed to be the saae as those 
at design condition. 



OHaPT^H III 


MTHQD OF SOLUTION 


III .A. FIIUTg DIFF^fidlNCi g ■^.UaTXOINB ; > 

In order to solve the radial equilibrium equation 
01 , 19 ^ and the continuity aquation (ll, 23 ), a f i nite -dif fe re n 
method is used. The fi nite -differs nee grid in the meridional 
plana is shown in Fig, 1, 

The first-order derivatives with respect to stream- 

8C|. ^ , 2 3 jj 

line radius, 'r^, like — , , -r — etc., which appear 

d r Cl r or 

in the radial equilibrium equation, are approximated by the 
backward difference. That is, at grid point (i,j), 




3 r ■ 


3-1 j 








(III, l) 


v^ith the help of the above equation the radial equi- 
librium 6quation5^ while integrated in radial direction, for 
a compressor is transformed into the following f i nite -dif f e re 
e quatio n. 
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For a turbine n is replased by 1 / n , . 

c t 

The equations (ll, 25 , II »29 and II« 30 / in the finite 
difference scheme are as followss 


i. . = H. 

3 - 


^.+tD(r- .G„ ^r. ) 

-1,J 0i_j 


1 ° % 1 
XjJ 


Hence equation (lll« 3 ) becomes 


r-T = H. + 0 ) (r. . G tan a. . - r. , . -q 

‘^i,0 i-1,j ^ 3.,j x^^^ 1,0 1-1,0 6 


i- 


tan g- . = (Gq r. j ) / G 


Us i ng e qu atio ns (TII# 4 ) and( HJ »5 ) in equation (III« 2 ) 
following quadratic equation in 

i,j 


(111. 2 ) 

(111. 3 ) 
(I11.4i 

) 

1 » j 

(111. 5 ) 

(111. 6 ) 

, the 


is obtained 
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A O +BG +d^=0 


(III. 7) 


where 
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(2 +\ i + . j + r (n, + — ~)i X. X. . 

4^0 , ,.2 "'> X»J X,J 
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{III.S 


3 = - 0 ) tan a-^ 


(III. 9) 


and C = .. 
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the integral forn is 

expressed in the follow! ng f ini te -diffe re noe form. 
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The siumerioal solution of tbs f i nits -dif f e re oce 
equations (HI. 7;* and (ill, with ths help of equations 
(lll,4j to (ill, 6) is described in the next section, 

III.3, SOLdflOH OF FUgfiS ^DIFF^H.! 

The calculating grid in the meridional plane is shown 
in Fig,1, Calculating stations are taken after rotor and 
stator blade rows. For calculating stations after rotor blade 
rows 0 ) is the angular speed of ths rotor and for calculatirg 
stations after stator blade rows » is zero. The exit flotv at 
each blade row is taken as the inlet flow to the succeeding 
blade row. The calculation proceeds in the following steps, 

1, At each station the annulus is first divided into a number 
of streamtubes of equal area. In this case l/£/^xe numbe r of stream 
tubes, where 'f^ is the fraction of mass flow passing through 
each streamtube. The streamlines, thus generated serve as 
temporary guides in the iteration process and are repl ajed by 
newly calculated streamlines af te r e ach overall iteration, 

2, The first three caloul ating statio ns are take n outside the 
first blade row. The flow field at two stations is calculated 
assuaiirg The slope and curvature of 

d( 

streamlines of these two statio ns are^ to be zero. From third 
station onwards the calculation proceeds as follows. 
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3. The propertiss of streanlines i*6, elope and curvature at 
all grid points are computed by fitting a tbres-point parabola. 


4 . For the first overall iteration, the axial velocity at 
each grid point is assumed to be equal to the axial velocity at 
the first two stations. 

5 . each grid point, the terms such as M, A , S , Gg , I-I, p and 

tan B are 4^ — calculated from assumed 0 and given 

V alues of tan Cl , 


6, The radial equilibrium equation is nov; integrated from 
streamline (j-l) to streamlins new value of C is obtained 

by solving the quadratic equation (III* 7 ;« 


7 , The nevj value of 0 is checked f or c onve rgs nee , If it 

is not converged to a specified tolerance, steps 5 and 6 are 
repeated us i the following value of axial velocity. 


G. 


X. 


= G. 


If J 


^k-l 

'x. . 


( r,k 

V 

' X. 


J 


k X 


( we ighte d i 


(III, 12 / 


where 

rs 1 axat io n 
ve lo city c 
error. For 


is the relaxation factor on axial velocity. This 
factor is necessary to avoid divergence of axial 
alculations* The value of is found by trial and 
the present case is taken as 0* I4 
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8 , After axial velocity is cslculated at each grid point 

in a particular calculating station, new values of streamline 
radii at that particular station are calculated from the 
continuity equation 

9 , The tip radius, thus calculated, should coincide with 

the given tip radius, at that particular station, 3ut, usually, 
a tip radius, woich is either small or large conpared with 
the given tip radius is calculated, new value of hub axial 
velocity is calculated as follows and used in the next iter- 
ation to calculate tip radius. 

For the second iteration, the hub axial velocity is 
taken as O .9 times of the hub axial velocity of first itera- 
tion. From third iteration onwards, the hub axial velocity 
is determined by linear interpolation of hub axial velocities 
of the previous two ite r atio nit. Thus, for third iteration, the 
hub axial velocity is given by intersection of ^ 

and for n • iteration, it is give n by intersection of 
P^ ^ and A3, as shown Fig, 2. 

10, Steps 5 to 9 are repeated with the new bub axial velocity 
until the calculated tip radrus agrees with the given tip 
radius within specified limit, 

11 , In this way, the solution marches down uptc 
oalculatir^ station, which is the- first station outside the 



26 


last blade row. Here M is the total nunibe r of calculatir^ 

stations, J^it the last two stations, the flow properties are 

th 

assumed to be the sans as those of the (m- 2 } station, 

12, For the next overall iteration, new streamline radii are 
obtained by adding a 'weighted fraction of the shift in radius, 
to the radii of the previous overall iteration. That is, 

r. . , V ~ ^ . V + (.r. . - r ) 

(weighted) 2 x,j i,j-l 

where is the relaxation factor on streamline radial shift, 

15 

^'s/'ilki ns on derives explicit expressions for the optimum rela- 
xation factor as function of grid aspect ratio, Mach number 
and the method for finding slopes and curvature of the meridio- 
nal streamlines for a simple two-dimensional flow. For the 
present study of axisymme trio flow, explicit expression for 
relaxation factor is not derived. It is take n to be 0,1 as 
suggested in reference 17, 

13 * Using the weighted streamline radii, the calculaticn from 
steps 3 to 12 is repeated* 

14, In suGcessive overall iterations, the calculated stream- 
line radii converge within a specified tolerance , Us u ally 10 
to 15 overall iterations are tjeeded for converge ncy on stream- 


line radii 
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15. “^‘fter calculating the velocities and other flow 
props rtias at e ach grid point, the mass averaged psrfor- 
aianos after each stage is calculated using the procedure of 
re fe re nee 9 , 

The mass»ave rage d temperature rise across a blade 
row is given by the following equation* 


2TT K 


0 


^i,tip T 

/ (-T 

^i,hub * 


I III. 13) 


where is, the inlet stagnation temperature and is tem'* 
pe rature-rise or energy exchange factor. This factor is 
necessary, in addition to the boundarj' 1 r blockage factor 
, to account for the wall boundary layer effect on tsapera.- 
ture - rise , 


Similarly, a correction factor for the mass*- 

P 

avs raged pressure ratio is applied and the mass -averaged 
pressure ratio is given by the following equation.' 

P. K p, Idll _ 

{■—y = -_£- {/ [(^) - 1 ] P dr^}+ 1 (III, 14) 

^ ^i.bub 


where Tp is the inlet stagnation pressure . 

From survey of limited experimental data, the values 
of i£ and & are taken to be 1. i n re f e re nee 1D 

P B ' 
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After niass*»averaged pressure ratio and temperature rise, the 
adiahatio efficiency and polytropic efficiency are determined stage- 
hy-stage, 

1=1 

(p./p ) ^ - 1 

Adiahatio Efficiency « — - — (ill, 15) 

^ (T./T ) • 1 

.In (Pi/Po)„ , .X 

Polytropic Efficiency n_ « i~^) ^ (ill, 16) 

^ ^ In (T./t ) 

In the following section details of numerical difficulties 
encountered in calculation of velocities and streamline properties 
and their elimination procedures are described. 

II. C, NUMERICAL DIFFICULTIES AM) THEIR ELIMimTION s 

The accuracy of any streamline curvature method depends on the 
accuracy of estimating streamline properties, i.e,, slope and 
curvature of meridional streamlines. Since, during successive 
iterations, the streamline may not be smooth, a suitable method 
of curve fitting should be employed. Especially the curvature 
is very sensitive to the smoothness of the streamline. In the 
present study, the spline method described in reference 29 
is first employed to find streamline slope and curvature. 

Even though calculated ' slope is reasonably accurate, . curvature 
calculated by this method is usually Very high, especially when 
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tbe atrsamiine ig not smooth. This high value of curvature 
gives difficulties in calculating axial velocity. Afterwards 
a second degree polynominal passing through three successive 
points is employed. This method gives fairly accur atsjs lope 
curvature, even when the streamline is not smooth. 


a nd 


Daneshyar and Shalaan had shown that under certain 
conditions tbe truncation errors in the numerical solution can 
become large and hence give rise to the violation of uniqueness 
conditions. The •'^um^rioal method may then give wrong answers 
to the physical problem. The conditions for existence and 
uniqueness may be violated when axial velocities are small 
(e.g, near staliy or when there ere regions of choked flow. In 
this case the program may give large value of tip radius. Hence 
in adjusting tip radiug, small velocities may occur at some 
streamlines. This may further give rise to negative velocities. 
Since negative velocities give difficulties in calcul sting other 
flow parameters, these velocities are replaced by some reason^l 
positive values , 

In determining axial velocity from radial equilibrium 
8 quatio n by the quadratic equatio n ( 111,7/, the discriminant 

' p' ' ' . 

( i.e , B -4AG/ may b eco me negative , He nee solutio n of r,e ,e . f ail 
From expression for A (III*8), it is seen that A is Bl^ays \ 
positive. Hence should be small or negative, so that the 
discriminant is positive, Tbe value of S mainly depends upon 
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absolute oi rcuoif e re nti al velocity. Hs ncs tbs discrisi nsbt 
be coaes negative « te n eith = r G or a is lar^e, 3o the value 
of O^is limited to a maximum value. This further eliminates • 
very high Mach numbers, High Mach numbers may give riseU? 

p 

numerical instability because the term l/( 1-M > appears in 
r,e,e, vs n whs n 0! is limited, the dis crimi nsnt may become 

JX, • 

negative. In such cases, the discriminant made equal to zero 
for small negative values. For high negative valnec c:? c:lc- 
criEiinant, the progras terminatsE present data anc .proeo-dc 
OQ to next data. 
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RESULTS MD DISCaSSICM 


The method described in previous chapters is used to predict 
the performance of the following turbomachines. 

26 

1. Twowstage Transonic .Axial-Tlow Fan 

27 

2. Five-stage Transonic Jixial-Elow Oonpressor 

28 

3. Two-stage aibsonic Axial-Flow Turbine . 

Before going into details of results and discussionj theoretical 
variation of perfonnance parameters with axial or meridional velocil5^ 
is described. 

2 

Meridional velocity 1+ X 17.1 

is approximately equal to for small values of 

Absolute Mach number is given by 

M . = / (l + X^ + tan^ a) TJ,2 

abs a 

and speed of sound’ 

= (1 + a)]l7.3 

Since is always smaller than H, 'a* is approximately constant 
X 

for small changes in 0^, and o at desigji point is same as that in 
experiments. Hence radial variation of essentially follows that 

of 0^ or 0^. 



Rslative Mach nxambcr is given by 


= (1 + H. tan^ <. . 22£-iSa_“ + iL^) 


If the Y/heel speed tor is of higher order than then 

discrepanoes between experimental and theoretical values of will 

be small otherwise, the radial variation of M . follows that of C . 

' rel X 

Relative fluid angle for a rotor is givai by 


6 m tan* (tan a — 

0 

X 


Hence g decreases as increases and converse is true. 


Total temperature ratio across a rotor blade row ’i^ and 
along streamline ' j ’ is given by 


T. ^ . H. ■ . 

1-1 »J 1-1 #3 


^i-1,3 


Since r. . and r. ^ . are approximately equal, 

1>3 1-1,3 


T. . mr. . ^ 

^ “)i J - (0^ tan n)^. ] 

i-1 ,3 1-1 r3 


Henoe total ten5)era1aire ratio is dependent on three parameters, 
avT velocity at leading and trailing edges of rotor and leading edge 
stagnation enthalpy. 33ae total pressure ratio variation essentially 
follows that of the total toaperatare j^tio. 
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In the following sections, overall performance and radial 
variation of performance parameters are presented for the three 
examples . Also consistency in the predicted and expearimental 
variation of various parameters is examined and wherever possible, 
explanations are offered. 

IV A. iPwo-Stage Iransonic Axial-Plow jPan (Beference 26 );- 

Ihis fan was desigaed for a tip speed of 1450 ft/sec, an overall 
pressure ratio of 2.8 and a corrected mass flow of 184.2 Ib/sec. Eie 
fan was designed without inlet guide vanes and air enters the first 
rotor axially. Eie inlet hub-to-top ratio was 0.4- and exit hub-to- 
top ratio was 0.68. Rotor and stator blade sections for both sta-ges 
of the fan are multiple-circular-arc (mCA) airfoils designed on 
conical surfaces which approximate stream surfaces of revolution. 

!Ehe computing time taken on IBM 370/165 computer for the 
execution of oatire performance is about 2.5 nonutes. Por a 
typical point near desigi condition, the computing time is 10 to 12 
seconds. 

IV. A. 1 Overall Performance The overall performance data for 

the two-stage fan are presented in figure 4. These data cover a 
range of corrected speeds from 50 to 100 percent desigi. The pressure 
ratios predicted by this method are usualdy higher than the experimental 
values. The efficiences are in good agreemmt with the experimental 
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values, whicii is to be expected, since the actual value of polytropic 
efficiency at each point is fed in as input. 

Miabatic efficiency at desiga point is 81 percent, which is 
almost equal to the experimental value of 80.8 percent. Peak 
efficiency at desigi speed is 84.4 and occurred at a pressure ratio 
3.29. Corresponding experimental values are 85 percent and 2.95 
respectively . Peak efficiency increased from 82 percent at 50 percent 
design speed to 84.4 percent at design speed. Experimental peak 
efficiency at 50 percent desiga speed is 82.2 percent. 

Overall performance of first stage is presented in Eigure 5. 
Efficiency at design operating condition is 81.8 percent which is 
3.6 percent less than the experimental value. Pressure ratios at 
desiga speed is in good agreement with eaperimaatal values. At other 
speeds the predicted values of pressure ratio are not in good agreement 
with the experimental values. Ihe predicted efficiencies are more 
scattered than the experimental values for the first stage as oompar«d 
to the whole two-stage fan. This may be due to the assumption of 
constant polytropic coefficient for all stages, 

I7.A.2 L'Spaawise Variation of Pertinent Parameters at Design Point . 

The predicted spanwise variation of pertinent parametei« at 
desigi point for two-stage fan is presented and compared with 
experimental variation in figure 6, 
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The meridional velocity profile at"Lr£.of rotor 1 predicted 

by this method is higher than the experimoital profile, Spanwise 

variation of M , , M _ and B are consistent with C variation, Eo 
abs rel m 

definite explanation is offered for the differences in profiles, 
but it may be due to assumption of one-dimensional flow and zero slope 
and curvature at the two stations before rotor 1 li.E. 

The predicted spanwise variation of at I.B. stator 1 is 
lower than experimental variation and at hub and tip regions the 
velocity profile tends upwards. This may be due to assumption of 
constant polytropic coefficient. This discrepancy can be checked 
by running the program with actual value of polytisopic coefficient 
at hub and tip. The differaices between spanwise variations of 
predicted and experimental absolute Mach number follows that of 
as expected. The I.E. ■ conditions of stator 1 are taken to be 
same as those at trailing edge (t»E. ) of rotor 1. Hence the total 
temperature and pressure ratio variations are lower than the 
experimental variations. 

At I.E. of Rotor 2, the a profile is lower than exp.erimental 
profile, ^anwise variation of other flow pioperties follows this 
trend. The total temperature and pressure ratios spanwise variations 
follows those of meridional velocities at l.B, of rotor 2 and 
stator 2* 
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IV. B. Bive-Stage Transoaic Axial-Plow Oompressor (Beference 27 ) 

This compressor was designed with all rotors operating with 
transonic relative inlet Mach numbers. The compressor was desigied 
for radially constant axial velocity and other parameters. Inlet 
hub-to-tip ratio was 0.5 and exit hub-to-tip ratio was 0.816. Ho 
inlet guide vanes are used. The stators were designed to turn the 
air axially. Both rotor and stator blade sections were double- 
circular arc (boa) profiles. 

The computing time taken on IBM 370/165 Computer for the 
execution of entire performance map is 4 minutes. Bor a typical 
point near design condition, the computing time is 20 seconds. 

IV.B.1 Overall Performance 

The overall performance data for the five-stage conqjressor 
are presented in figure 7. These data cover a range of corrected 
speeds from 70 to tOO percent design. The total or stagnation 
pressure ratios predicted by the present method are 5 to 10 percent 
lower than the experimaatal values, Bor this compressor a constant 
value of poly tropic coefficient is used for all operating points. 
Hence this results in a flat curve for the adiabatic efficiency. As 
expected, the efficiencies are in close agreemoat. The predicted 
efficiencies are tisually 2 to 3 points lower than the experimental 


values . 
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Oversll perf oimanc© of individual stages is not piesmted for 
this compressor. 

IV«B»2 Spanwise Variation of Pertinent Parameters at Sesi^ Joint : 

The spanwise variation of performance parameters at desiga 
point for this compressors are presented along with experimental profiles 
in figure 8, The spanwise variation of other parameters follov^s that of 
^m expected and no inconsistaacy is observed in either 

predicted or experimental profiles. But for later stages the variaticm 
in other parameters does not follow that of C , since at these stages, 
these parameters are affected by the discrepancies between exp erimentals 
and. predicted performance of first stages. This can be explained as 
follows 

TotaJL enthalpy H = H + «<j((rO tan a). . - (rC tan a). . .] 

AyJ yj * » yJ 

... iv.s 

Hence the discrepan(^ in values of first stages affect the 
value of H at later stages. Since H enters in calculation of other 
parameters, the differences between predicted and ex^serimaital values 
are pronounced in later stages. 

IV, C. Two-Stage Subsonic Axial^Mow Turbine (Beference 28) :- 

The turbine, having an overall loading of 3.25 for the two 
stages, was fitted with aerofoil camber— line typ6 blaaing. The base 
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profile shape was placed around a parabolic cambK^ line with the position 
of maximum camber at 40 per cent of the chord from the leading edge. 

The annulus was flared, with hub/to-tip ratio falling more or less 
linearly from 0.8 at inlet to the first-stage stator to 0.66 at exit 
from the second-stage rotor. 


The computing time taken on IBM 370/165 Computer for the execution, 
of entire performance map is 2 minutes, for a typic^'l point near desiga 
condition, the computing time is 10 seconds. 

IV.G.1 Overall Performance ;« 


aiie overall perfo^rmance data for two-stage turbine is 
presented in figure 9« Ihese data are presented at four speeds, 50 
percent (i.e, = 152), 75 percent (i.e. = 194), 87*5 percent 

(i,e. » 230 ) and 100 percent (i.e, = 263) desigi speed. 

The adiabatic efficiencies are in good agreement with the experimental 
values, as to be expected. The total pressure ratios predicted by 
the present method are about 10 percent below the experamaatal values. 
Horo the total pressure ratio is the ratio of total pressure at inlet 
to the turbomachine, to total pressure at exit of the turbomachine. ^ 



4-^;) 4 0 -fiPiire 10. Experimental data 
These data are presented in iigu 

. ^ vo- c+ato-r 1. rotoib 1 and stator 2 

are not available for compariscai. 

4 .S n-rroximately ccaietant along the 
the L.B. meridional velocity profile is ai^prox ^ 
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spaa. Other parameters follow this tread. Star rotor 2, the 

20 

variatioa at tip is marked. Erost had obseirved same trend at 
87 .5 percent design speed. 

Prom the ahore discussions, it is observed that the spanwise 
variation of various parameters closely follows that of meridional velo- 
city. Once the discrepancies between predicted and exp. erimental 
values are explained, differ<aices in other parameters can easily 
explained. The reason for the discrepancies in may be due to 
assumption of radially constant polytropic coefficient among other 
factors. This can be confirmed, if the program is run v/ith variable 

poly tropic coefficient. 
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CONdDSIOMS 

A computer program for predicting the perfoimaice of turbo- 
machines has been written. This has been done by taking calculating 

stations after blade rows and solving for axial velocity and stream- 
line radius. 

The following turbomachines have been selected as test cases 
for the program; a two— stage transcnio axial— flow fan, a five— stage 
transonic axial— flow compressor and a two— stage subsonic axial— flow 
turbine. 

!I!he program predictions have been compared with experimental 
data, except for two-stage tu 3 *ine. 

The results obtained by ■fee use of the present program are 
useful as a first approximation to determine the performance of 
turbomachines. The computed overall pressure ratios agree to within 
5 percent of the experimental values in almost all cases tested, itor 
compressors, the approach of stall is correctly predicted by the 
computed performance map, i^ich is a useful feature. 

The comparisons with the experiment show that this method gives 
a reasonable estimate of spanwise distribution of flow properties. An 
improvement on the measure of agreement demands a better model for the 
representation and distribution of losses. 

The program was not run over a sufficiently wide range of parame 
ters such as H, ^ deteimine effect of these parameters. 
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SUGGS STIOIg FOR FWTHBR ^ORKt 

In the Present study, th s po Ly tropic coefficient is 
assumed t o be constant at all grid points , This means that 
losses are equal at all streamlines. This assumption is invalid, 
since in hub and tip regions losses are more than in mid-span 
region of tne bl aie , He nee , v ari ab Is poly t ro pic coefficient 

may be used to obtain more realistic performance data, 

3 

Swan developed correlation of loss parameter with 

equivalent diffusion factor D for D,0,A, rotor blaSe elements 

e q 

as a function of the percent of the span. These correlations 
may be extended to other types of blede and used in the present 
study. Then polytropic coefficient need not be fed in as input. 

For transonic compressors, shock losses may be significant. 
In the present study, these losses are neglected. In further 
work, shook losses may be taken into account assuming the normal 
shock model developed in reference 30 , 

lau and Lewis'"^ has recently demonstrated that sweep 

and dihedral of the blades can very significantly affect cascade 

32 

characte risij2'(^0 - . Smith and Yeh h ave airs ady de vslope Q a 
mode 1 f or i ncludi og tbs above effects. The present program 
could be improved by incorporating suitable corrections. 
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la the present study, the relative fluid angles at off- 
design conditions are detexinined using correlation of deviation 
angles with equivalent diffusion factor. lo obtain more realistic 
cascade characteristics ^ analytical cascade model should be developer 
and included in the present program. 

3h all mathenatical models of tufbcanachine flow field, the 
formulation starts with inviscid fluid assumptions, but losses are 

introduced by means of a poly tropic coefficiaat or pressure loss 

5 

coefficient. Sbiith had reserraticns about this approach. 

22 

Horloofc had takm into account of this fact and formulated the i 

21 i 

prohlem starting with viscous flow equatims. Bosman and Marsh i 

had included this correct icaa in their foimulation for matrix 
throng flow method. She present foimulaticaa may be modified in 
those lines. Jlnally any theoretical method should be compared TO.th 
other theoretical aethods. Saaee the reazlts of this program may 
be compared with those of references 19,20 and 21. 
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